Quantum discord and entanglement in Heisenberg XXZ spin chain after quenches^ 



Ren Jie, 1 WU Yin-Zhong, 1 and ZHU Shi-Qun 2 

1 Department of Physics and Jiangsu Laboratory of Advanced Functional Material, 
Changshu Institute of Technology, Changshu 215500, P. R. China 
2 School of Physical Science and Technology, Soochow University, Suzhou, Jiangsu 215006, P. R. China 

Using the adaptive time-dependent density-matrix renormalization group method, the dynamics 
of entanglement and quantum discord of a one-dimensional spin- 1/2 XXZ chain is studied when 
anisotropic interaction quenches are applied at different temperatures. The dynamics of the quantum 
discord and pairwise entanglement between the nearest qubits shows that the entanglement and 
quantum discord will first oscillate and then approach to a constant value. The quantum discord 
can be used to predict the quantum phase transition, while the entanglement cannot. 
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Entanglement is an important resource in quantum in- 
formation processing, such as quantum computation and 
quantum teleportation [TJ [5] . Many researches show that 
entanglement exists naturally in spin systems at zero 
temperature. By using the quantum theory, entangle- 
ment can be applied to detect the quantum phase tran- 
sition points [3H5] , which imply quantum fluctuation due 
to the change of a parameter in the Hamiltonian at zero 
temperature. Recently, a different and significant mea- 
surement, called quantum discord, was introduced [TJ. It 
can detect these quantum correlations presented in cer- 
tain separable mixed states, which cannot be captured by 
the entanglement of formation. The discord can also be 
used to predict quantum phase transition points at zero 
temperature [81411] and at finite temperatures [12] . It is 
reported that the measurement of quantum discord can 
be realized in NMR [T3] and optical experiments [13] . 

Recently, the progress in manipulating cold-atom sys- 
tems provides almost ideal realizations of strongly cor- 
related many-particle systems in experiments [15] . The 
non-equilibrium physics in quantum many-body spin sys- 
tems are attracting much interests of experiments due to 
the possible application to spintronics, quantum informa- 
tion processing, etc. Meanwhile, physicists try to under- 
stand the general aspects of non-equilibrium dynamics 
in a comparably simple system. One of the methods is 
the quench where the system's Hamiltonian is changed 
suddenly. It corresponds to preparing the system in the 
ground state of the Hamiltonian Hi at the initial time, 
and then the Hamiltonian is changed instantaneously to 
Hp after time t > 0. The quench dynamics of entangle- 
ment and quantum discord in a spin- 1/2 transverse XY 
chain was studied [161 117] . Non-equilibrium dynamics of 
the XXZ model with quenches has been a very active area 
of research [TBrl2"2"] , since such a global interaction quench 
can be actually realized for atoms in optical lattices [23] . 
It is interesting to investigate the dynamics of entangle- 
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ment and quantum discord in Heisenberg spin- 1/2 XXZ 
chain following quenches. 

The Hamiltonian of an opened N spin-1/2 chain with 
the nearest neighbor interaction is given by 

JV-l 

H=J2 J(S*S* +1 + SfSf +1 + AS*S* +1 ), (1) 

8=1 

where Sf(a = x,y,z) are spin operators on the i-th 
site, and the parameter N is the number of spins in 
the chain. The parameter J refers to the ratio of the 
nearest-neighbor coupling. For simplicity, J = 1 is con- 
sidered. The parameter A denotes the anisotropic inter- 
action. The anisotropic interaction A > 1 is called Neel 
antiferromagnetic phase, and < A < 1 is called critical 
XY phase. In this paper, we assume that the anisotropic 
interaction is time dependent. The initial anisotropic in- 
teraction is labeled as Aj. After t > 0, the anisotropic 
interaction A/ suddenly quenches to Ap- 

In the Letter, the concurrence is chosen as a measure 
of the pairwise entanglement [23J- The concurrence E is 
defined as 

E = max{Ai - A 2 - A 3 - A 4 , 0}, (2) 

where the quantities Aj(i = 1,2,3,4) are the square 
roots of the eigenvalues of the operator g = p\i(?X ® 
a\ )p\ 2 { (J \ ® a V) an d m descending order. The case of 
E = 1 corresponds to the maximum entanglement be- 
tween the two qubits, while E = means that there is 
no entanglement between the two qubits. 

The quantum discord which measures the amount of 
quantumness in the state was introduced by [TJ, 

Q(p)=I(p)-C(p), (3) 

where I represents the total correlation and C is the 
classical information. For an arbitrary bipartite state 
Pab, the total correlations, which are expressed by quan- 
tum mutual information I(pab) in Eq. ([3]), is given 
by I{pab) = s(p A ) + s{p B ) - s(pab), and s(p) = 
—Tr{p log 2 p), where pa{b) is the reduced density ma- 
trix obtained from pab by taking the partial trace over 
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the state space of subsystem B(A). Quantum condi- 
tional entropy can be defined as s(p\{Bk}) — J2kPk s (Pk)- 
The final state pk of the composite system, which is the 
generalization of the classical conditional probability, is 
given by pk = ^-{1 ® Bk)p{I ® Bk), with the proba- 
bility pk — tr(j eg) Bk)p(I ® Bk) where / is the iden- 
tity operator of the subsystem A. The measurement 
based quantum mutual information takes the form of 
I(p\{Bk}) — s(p A ) — s(p\{Bk})- The classical informa- 
tion is given by 



C( P ) 



J {B k } 



I(p\{B k }). 



(4) 



In Eq. ([3]), Q = means that the measurement has ex- 
tracted all the information about the correlation between 
A and B and is classical correlation. While Q ^ means 
that the information cannot be extracted by local mea- 
surement and the subsystem A gets disturbed in the pro- 
cess, a phenomena not usually expected in classical in- 
formation theory. 

Because of U(l) invariance [H, J2toi $ z ] — 0, the ex- 
pectation values (S z ) = in the ground states of Eq.([T]), 
and the evolution of the reduced density matrix can 
be described by p(t) = \[I AtB + ^ (a k A a k B (t))a k A a k B }. 

k—x,y,z 

The entanglement in the reduced density matrix p(t) can 
be express by [251 [25] 



£(p(t)) = imax{0, £ 

k—x,y,z 



(5) 



and the quantum discord is given by |27] 



Q(p(t)) = [5(1 

g(l - d x 4 
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with dk(k 
m.ax{\d x \, \d y \, 



(6) 
D = 



= x,y,z) = (<J k A a k B {t)} 
\d z \}, and g(x) = xlog 2 {x). 
The dynamics of entanglement and quantum discord 
of a spin-1/2 antiferromagnetic Heisenberg chain is an- 
alyzed at zero temperature. Two kinds of initial states 
are selected. The first initial state is chosen in XY phase 
region, and the second kind of initial states is in Neel 
phase region. We can use exact diagonalization to ob- 
tain the ground state for a small system, and then ap- 
ply size scaling behavior. If the system size is large, the 
adaptive time-dependent density matrix renormalization 
group (t-DMRG) method can be applied with a second- 
order Trotter expansion of the Hamiltonian [28 , 29J. In 
our simulation, a Trotter slicing dt = 0.1 and the t- 
DMRG codes with double precision are performed with 
a truncated Hilbert space of m = 500. We calculated 
a chain of L — 60 sites with time t < 10 and the rela- 
tive error is below 10 -4 . In the simulation time region, 
our numerical results are reliable, which have been esti- 
mated by keeping trace of the discard weight, and the 
entanglement entropy increases linearly in time [191 130] . 
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FIG. 1. (Color online) The entanglement (a) and quantum 
discord (b) are plotted as a function of time for different Af, 
when A/ = and the temperature of system is zero. The 
entanglement (c) and quantum discord (d) are plotted as a 
function of Af at time t = 10. Here and in the following all 
quantities are dimensionless and the dots are t-DMRG data, 
the lines are fit lines. 



Firstly, the initial state of Eq. ([!]) is in the XT phase, 
then the anisotropic interaction quench is applied. The 
entanglement and quantum discord between the two cen- 
tral qubits are plotted as a function of time for different 
Ap in Figs. [TJa) andJIJb) respectively for A/ = 0. In Fig. 
[lja), when Ap — 0.5,1.0, the entanglement increases 
firstly with the time increases, then reaches a maximal 
value, and finally decreases to a constant value. When 
Ap = 2.0, 3.0, the entanglement will first decreases with 
the time t increases, then it drops to a minimal value, and 
finally reaches to a flat with small oscillations. In Fig. 
[TJb), the dynamics of quantum discord for Ap = 0.5, 2.0 
is similar to that of entanglement. For Ap — 1.0, quan- 
tum discord has more oscillations in the beginning. To 
investigate the effects of the anisotropic interaction, the 
entanglement and discord are plotted as a function of 
Aj? in Figs. [ljc) and lid) at the fixed time t = 10 
respectively. In Fig. fl^cj, the entanglement decreases 
with Ap increases. While in Fig. |TJd) , the discord in- 
creases slightly when Ap increases, and reaches the max- 
imal value 0.283 when Ap = 1.0, and then the discord 
decreases wihen the interaction A^ increases further. We 
also investigate the cases of Aj = 0.25,0.5, the discord 
also reaches the maximal value when Ap — 1.0 at the 
fixed time. 

Secondly, the initial state of Eq. ([I]) is in the Neel 
phase, then the anisotropic interaction quench is per- 
formed. The entanglement and quantum discord between 
the two central qubits are plotted as a function of time 
for different A^ in Figs. [2]ja) and [2jb) respectively for 
A j = 3. In Fig. [2ja), when Ap = 0.5, the entangle- 
ment decreases with the time increases, then tends to a 
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FIG. 2. (Color online) The entanglement (a) and quantum 
discord (b) is plotted as a function of time for different Af, 
when Ai = 3.0 and the system temperature is zero. The 
entanglement (c) and quantum discord (d) is plotted as a 
function of Af at the time t = 10. 



FIG. 3. (Color online) The entanglement (a) and quantum 
discord (b) are plotted as a function of time for different Af, 
when Aj = and the system temperature is T = 0.5. The 
entanglement (c) and quantum discord (d) are plotted as a 
function of Af at the fixed time t = 5. 



constant value. For Ap = 1.0,1.5,2.0, the entanglement 
increases with the time increases, then oscillates to reach 
a constant value. In Fig. [5Jb), the discord increases with 
the time increases, then oscillates and approaches to a 
constant value. To find the effects of the anisotropic in- 
teraction, the entanglement and discord are plotted as a 
function of Ap in Figs. [2jc) and[5Jd) respectively at a 
fixed time t = 10. The entanglement increases when the 
interaction Af increases. It reaches the maximal value 
0.327 at A P = 1.7, and decreases with A F increases fur- 
ther. In Fig. [2^d), the discord increases linearly with 
the interaction increases, and reaches the maximal value 
0.271 at Af = 1.0, then decreases with Af increases. 
We also study the cases of A/ = 2.0, 2.5, the discord also 
reaches the maximal value at A^ = 1.0 at the same time. 

The dynamics of entanglement and quantum discord 
of a spin-1/2 antiferromagnetic Heisenberg chain is in- 
vestigated at finite temperature. It is also started from 
two kinds of initial states as before. It is known that 
the density matrix for a system in equilibrium at a fi- 
nite temperature T can be expressed by p — e z , with 
j3 = 1/kT and Z — Tr(e~@ H ) is the partition function. 
The Boltzmann constant k is set to unity. It is well known 
that the finite-temperature density matrix renormaliza- 
tion can be calculated by using an enlarged Hilbert space 
with imaginary time evolution [31] . , and this method can 
be extended to study the evolution of the system at finite 
temperature [32]. A simply method, the time-reversed 
physical Hamiltonian is applied to reduce the errors |33j . 
Here, we adopt this method to study the dynamics of 
entanglement and discord. In our simulation, a Trotter 
slicing dt = 0.1,(1/3 — 0.05 and the finite temperature t- 
DMRG codes with double precision are performed with 
a truncated Hilbert space of m = 300. We calculated a 



chain with L = 30 sites within t < 5, and a typically dis- 
carded weight can be kept below Sw — 10~ 6 with relative 
error below 10~ 4 . 

Firstly, the initial anisotropic interaction is selected as 
A j = and the temperature is chosen as T = 0.5. The 
entanglement and quantum discord between the two cen- 
tral qubits are plotted as a function of time for different 
Af in Figs. [3]ja) and[3^b) respectively. The phenomenon 
are similar to the case at zero temperature. However, 
when Ap = 2.0,3.0, the entanglement decreases to zero 
rapidly with time increases. If time is fixed at t = 5, 
the entanglement and discord are plotted as a function 
of the different Af in Figs. [3^c) and|3jd) respectively. 
In Fig. |3jc), the entanglement decreases very fast and 
reaches zero at about Af = 1.5 when Ap increases. In 
Fig. ^d), the discord decreases very slowly for Af < 1.0. 
When Af > 1.0, it drops quickly to very small value as 
Af increases. Though there is no peak at Af = 1.0, the 
sudden change of the slope at Af = 1.0 can be a hint 
about the quantum phase transition. 

Secondly, the initial anisotropic interaction is A/ = 3.0 
and the temperature is T = 0.5. The entanglement and 
quantum discord are plotted as a function of time for 
different Af in Figs. [4^a) andQb) respectively. They 
are similar to the cases of the quench at zero temper- 
ature. The entanglement and discord are plotted as a 
function of Ap in Figs. |3Jc) and[4]^d) respectively when 
the time is fixed. The phenomenon are similar to the 
case of the quench at zero temperature. The entangle- 
ment have a peak located at Af = 1.7 for T = 0.5, 1.0, 
and disappears at T = 2. The quantum discord have a 
peak located at Af = 1.0 and the peak value decreases 
with the temperature increases. 

At zero temperature, the dynamics of the entangle- 
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FIG. 4. (Color online) The entanglement (a) and quantum 
discord (b) are plotted as a function of time for different Af, 
at temperature T = 0.5 for Aj — 3.0. The entanglement (c) 
and quantum discord (d) are plotted as a function of Af for 
different temperature at t = 5. 

ment shows whether there is a peak in the entanglement 
with fixed time depends on the initial anisotropic inter- 
action. While for the dynamics of quantum discord, it 
reaches a maximal value when the final anisotropic inter- 
action is one, irrespective of the initial anisotropic inter- 
action. Therefore this behavior can be employed to pre- 
dict the quantum phase transition from the critical XY 



phase to the Neel antiferromagnetic phase [IS]. The en- 
tanglement fails to predict the quantum phase transition 
point while the quantum discord could detect it. At finite 
temperatures, when the initial anisotropic interaction is 
smaller than one, the quantum discord only provide a 
hint of the corresponding quantum phase transition. It 
is due to the very blunt of the peak of the quench when 
the temperature is zero. When the temperature increases 
a little bit, the peak disappears but a sudden change of 
the slope in quantum discord still exists. When the ini- 
tial anisotropic interaction is much larger than one, the 
quantum discord can predict the quantum phase transi- 
tion, because the peak of the same quench is very sharp 
at both zero and finite temperatures. 

To summarize, the adaptive time-dependent density- 
matrix renormalization-group method is used to investi- 
gate the dynamics of the entanglement and quantum dis- 
cord in one-dimensional quantum spin systems when the 
global anisotropic quenches are performed. The entan- 
glement and discord oscillate to reach a constant value, 
which depends on the initial and final anisotropic interac- 
tions. At zero temperature, whether there is a peak in the 
entanglement depends on initial anisotropic interaction. 
While for quantum discord, it reaches a maximal value 
when the final anisotropic interaction is one, no matter 
how the initial anisotropic interaction is. This behav- 
ior can be used to predict the quantum phase transition 
points. At finite temperatures, the discord can also be 
used to predict the quantum phase transition only when 
the initial anisotropic interaction is greater than one. 
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